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THE CALCULATED RESPONSE OF A THREE-CRYSTAL 
PAIR SPECTROMETER 

By Clemans A. Powell, Jr. 

Langley Research Center 

SUMMARY 

An investigation has been made of the use of a three -crystal sodium iodide pair 
spectrometer for the quantitative measurement of photon spectra over an energy range 
from 1.5 MeV to 10.0 MeV. The investigation included developing a computer analysis of 
the response and efficiency of the detector system and using an iterative unfolding method 
to analyze pulse-height distributions obtained with the system. The efficiency analysis 
considers the specific geometry of the three crystals and the attenuation of photons in the 
housings and reflectors of the assemblies. Simplifying assumptions used are as follows: 
(1) The pair-production processes of the incident photons occur on the central axis of the 
center crystal and (2) the resulting electron-positron pairs are brought to rest at their 
point of creation. The response of the spectrometer to monoenergetic photon sources 
was determined by using the efficiency calculations and the experimentally measured 
response to a known source. The iterative unfolding method uses these calculated 
responses to obtain experimental photon spectra from pulse-height distributions. 

The predicted efficiencies agree well with experimental efficiencies measured from 
an americium-beryllium source. The unfolding procedure has been tested by using sim- 
ulated spectra and has proved to yield errors much smaller than the errors introduced by 
either the predicted efficiencies or responses.. 

INTRODUCTION 

For many years the sodium iodide (Nal) scintillation detector has been the basis of 
various forms of spectrometers for measurements of gamma rays, characteristic X-rays, 
and bremsstrahlung. Among these forms are simple, single-crystal detectors; cylindri- 
cal crystal detectors with anticoincidence annuli; two crystals in coincidence; and three 
crystals in coincidence. All the latter forms are designed to overcome one of the major 
shortcomings of the simple, single-crystal detectors. Because of Compton scattering and 
losses of pair-production annihilation quanta, the response in the form of pulse-height 
distributions from a single-crystal spectrometer is quite complicated, especially when 
measuring energies above 2 MeV. However, the three-crystal spectrometer is capable of 


distinguishing between pair -production processes and the other forms of interaction of 
photons with the scintillation material and yields pulse-height distributions with single 
peaks for monoenergetic photons. Some previous uses of such a spectrometer are 
reported in references 1, 2, and 3. 

Before any type of spectrometer can be used in a quantitative experiment, the 
response of the spectrometer must be known over the photon- energy range of interest. 
Because there are only a small number of readily available sources of known photon 
energy above 1.5 MeV (the lower limit of the spectrometer's usefulness), a means of the- 
oretically predicting the response of the spectrometer is necessary. In addition, some 
means of unfolding or analyzing the experimental pulse-height distributions is necessary 
to yield meaningful experimental spectra. The present effort was made to furnish the 
needed information. 

Theoretical predictions have been made of the efficiency of a three-crystal pair 
spectrometer and its response to a range of monoenergetic photons, and an iterative 
method has been found to unfold the experimental pulse-height distributions. Not all 
entering photons are sensed by the spectrometer, and the efficiency analysis evaluates 
the effect of the geometry of the three crystals and the attenuation losses of photons in 
the reflectors and housings of the detector assemblies. The response of the spectrom- 
eter to a range of monoenergetic photons has been determined from the calculated effi- 
ciencies and the shape of an experimental spectrum of a known monoenergetic source. 
The unfolding method is based on an iterative procedure developed in reference 4. 

SYMBOLS 

B response matrix for spectrometer 

B element of B 

C matrix defined by C = BB 

D diagonal matrix; N = D S 

D element of D 

d distance along £-axis from point of annihilation-quanta creation to origin of 

coordinate system, cm 

E photon energy, MeV 


2 


F 

H 

I 

i,j,k 

L 

M 

M 

N 

N 

m,n,p 

P 

q 

R 


r 


S 

S 

T 

t 


u 


peak-to-total interaction ratio (photofraction) for side detectors 
height of side crystals, cm 

indentation of center detector into well bounded by side detectors, cm 

increment numbers for t, 6, and cp, respectively 

length of center detector, cm 

matrix of pulse-height distribution 

element of M 

matrix of photon-energy spectrum 

element of N 

indices of matrix elements 

probabilities of interactions making up resultant efficiency of spectrometer 

distance from point of annihilation-quanta creation to outer surface of side 
detectors, cm 

radius of detector, cm 

distance from point of pair production on Z-axis to surface of center crystal, 
cm 

matrix defined by S = BM 
element of S 

sum of squares of difference in Sp and its iterative approximation 
distance of point of pair production from origin along Z-axis, cm 
maximum value of m 
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V 

X,Y,Z 


X',Y' 


x 

X' 

y 

y’ 

a 

Aft 

A0,A<p 

At 

0 ’ 

P 

P 


maximum value of n 

coordinate system with origin at center of front face of center detector 

floating coordinate system with origin at point of annihilation-quanta 
production 

detector light shield thickness, cm 
annihilation-quanta path length through light shield, cm 
reflector material thickness, cm 

annihilation-quanta path length through reflector material, cm 

total number of increments 

incremental solid angle 

incremental polar angles 

increment of center detector thickness, cm 

coordinate system with origin at center of end face of side detector 

polar angles from point of annihilation-quanta production, rad 

polar angle measured from £-axis to a point on surface of side detectors, rad 

mass attenuation coefficient, cm _ l 

radius of hemicylindrical groove in face of side detector, cm 

distance from origin of coordinate system to a point on surface of side 
detectors, cm 


Subscripts: 
a aluminum 
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c center detector 

e total attenuation for Nal 

f stainless steel 

g pair-production attenuation for Nal 

i,j,k indices for increments of t, 9, and cp, respectively 

m,n,p indices of matrix elements 

r reflector material AI 2 O 3 

s side detector 

t thickness of center detector 

u maximum value of m 

v maximum value of n 

9,(p polar angles 

1,2 identification numbers of annihilation quanta 

Superscripts: 

q iteration number 

1 first iterative approximation 

A tilde (~) over the symbol for a matrix indicates the transpose of that matrix. 

SPECTROMETER OPERATION 

A sketch of the detector described in the present report is shown in figure 1. The 
spectrometer consists of a center Nal scintillation detector, 4.12 cm in diameter and 
10.16 cm long, which is coupled to a photomultiplier tube inside a magnetic- shielded 
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housing. In addition, there are two Nal annihilation detectors. Each detector is 
12.70 cm in diameter and 7.62 cm in height and has a photomultiplier tube and a magnetic 
shield with a hemicylindrical slot of 2.54-cm radius machined across the end face. 

During operation, the two side detectors are placed face to face, and the center detector 
is placed in the cylindrical well formed by the grooves in the side detectors. 

Photons incident on the center detector interact with the scintillation material (Nal 
activated with thallium) by three distinct processes: photoelectric effect, Compton scat- 
tering, and, if greater than 1.022 MeV, pair production. It is assumed that if a pair- 
production process occurs, the resulting electron-positron pair deposits its energy in the 
detector and is brought to rest. The positron combines with an electron in the detector, 
and the pair promptly annihilates, emitting two photons of equal energy (0.511 MeV) in 
opposite directions. If these photons escape being captured by the center detector and 
are captured by the side detectors, the events are recorded. Figure 2 shows a block 
diagram of the electronics used to insure that the necessary coincidences are met and to 
analyze the resultant spectra. The light output from the scintillation detectors is pro- 
portional to the energy deposited in each detector; consequently, if a pair -production pro- 
cess occurs in the center detector, the light output and electrical pulse from the center 
detector are proportional to E - 1.022, where E is the incident photon energy in MeV. 

In order to insure that the pulses recorded in the center detector were produced by a 
pair -production process, pulses from the center detector were recorded only when 
0.511-MeV pulses were registered in each of the side detectors. The single-channel 
analyzers, coincidence unit, and linear gate unit insure that only the proper pulses are 
recorded. 

The greatest benefit of a three-crystal spectrometer is the simplicity of the spec- 
trum produced by monoenergetic photons. Figure 3 shows the spectra produced in both a 
7.62-cm by 7.62-cm single-crystal Nal detector and the present three-crystal spectrom- 
eter by the 4.43-MeV gamma ray from the 9B e ( ajn y)12Q reaction. The single -crystal 
spectrum exhibits a complicated structure with three distinct peaks, corresponding to the 
full-energy peak and two escape peaks resulting from the escape of either or both of the 
pair -production annihilation quanta. There are a number of counts recorded in the lower 
channels because of losses from Compton scattered photons. The three-crystal spec- 
trometer spectrum, on the other hand, has a much simpler shape and is essentially a 
Gaussian peak plus a low-energy tail. This low-energy tail is the result of bremsstrah- 
lung losses and electron escape from the center detector. 

The primary disadvantages of a three -crystal pair spectrometer are low efficiency 
and because of the minimum energy necessary for the creation of the electron-positron 
pairs, a relatively high lower limit of photon-energy sensitivity. The rest mass energy of 
an electron-positron pair is 1.022 MeV, but the useful lower limit of photon-energy 
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sensitivity is approximately 1.5 MeV for . such a spectrometer because the pair- 
production cross section is extremely small at the lower energies. In order to obtain 
quantitative results from the use of a spectrometer as described in this report, whether 
it is used for measuring the spectra of discrete energy photons or continuous spectra as 
in the case of bremsstrahlung, it is necessary to know the efficiency and response of the 
spectrometer to monoenergetic photons over the energy range of interest. In addition, 
for continuous spectra or for spectra of many monoenergetic photons, a means of ana- 
lyzing the experimentally determined pulse -height distributions must be found. 


CALCULATION OF SPECTROMETER EFFICIENCY 


The efficiency of a three-crystal pair spectrometer is dependent on the specific 
geometry of the three detectors. Figure 4 gives the basic geometry of the detectors and 
the coordinate system on which the following discussions depend. This coordinate sys- 
tem has the origin at the center of the front face of the center detector and the Z-axis 
along the central axis of the center detector. 

In computing the efficiency of the system, the following processes are considered in 
turn, and the resulting efficiency is the product of the probability of the occurrence of 
each of them: 

(1) Passage of photons through the light shield and reflector of the center detector 

(2) Pair production at various points in the center detector 

(3) Capture of the annihilation quanta in the center detector and its light shield and 

reflector 


(4) Capture of the annihilation quanta in the light shields and reflectors of the side 

detectors 

(5) Capture of the annihilation quanta in the side detectors 

The resulting efficiency of the system will be shown to be of the following form: 


p[e] = P[entrance]^ P^[pair production]^ ^Pj k [AS2] p j ^ [escape] Pj k [pair capture]j 


where P[entrance] is the probability that a photon will pass through the light shield and 
reflector of the center detector, P^jpair production] is the probability that a pair- 
production event will occur at some position in the center crystal as determined by the 
indexing integer i, Pj^I^fl] is the probability that the annihilation quanta will be 
emitted in directions determined by the indexing integers j and k, Pj ^[escape] is 
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the probability that the annihilation quanta will escape being captured in the center crys- 
tal, and Pj } k[pair capture] is the probability that the pairs will be captured in the side 
detectors. 


Processes in Center Detector 

Gamma rays or photons impinging on the center detector from the negative 
Z -direction will be attenuated by the wall and reflector material surrounding the center 
Nal crystal. The material (not shown in fig. 4) used on the end window of the present 
detector is 0.04 cm of aluminum, and the reflector is 0.28 cm of powdered aluminum 
oxide ^A.l 2 C> 3 j . The probability of a photon entering the center Nal crystal is 

P [entrance] = exp 

where p a and p r are the attenuation coefficients for aluminum and the reflector 
material, respectively; and x c and y c are the respective wall thickness and reflector 
thickness of the center detector. 



All photons which interact with the center detector are assumed to enter and inter- 
act along the central Z-axis. This is a reasonable assumption, since, under operating 
conditions, a collimator is used to confine the incident beam to an area much smaller 
than the cross-sectional area of the detector. The point of interaction is the intersection 
of X' and Y', a distance t from the intersection of X and Y. If the displacement 
t is incremented into a number of thicknesses At, then 

t = i At (i = 1,2,3, ... a,.) (2) 


The probability of having a pair-production process in the ith thickness increment is 
given by 


Pj[pair production] = exp 


-p e (i ~ l)At 


1 - exp 


■p g 


At 


( 3 ) 


where p e and p g are the respective total and pair -production interaction (or attenua- 
tion) coefficients for sodium iodide. Pairs produced by Compton scattered photons are 
neglected. 

It is further assumed that the electron-positron pairs are stopped within the ith 
increment of their creation and that annihilation occurs in the same increment along the 
Z-axis. By assuming that the velocity vectors of the annihilation quanta created are uni- 
formly distributed, the probability of one of the gammas of a pair being emitted into an 
incremental solid angle AS2 with direction 9,(p and the other with direction 6 + tt, 
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cp + 7T is 


Pj jj[ao} ~ sin(j A9)A6 A cp 


(4) 


where 


0 = i Ae 


(p = k Acp - | 


i = 0,1,2, . 
\ 

k = 0,1,2, . 


.eta; A9 = — 
8 « 0j 


•v Acp = 2 cd 


Some of the quanta emitted from the ith thickness increment will be captured in the cen- 
ter detector. The probability that both annihilation quanta will escape from the center 
detector and its light shield and reflector is 


J*[§scape) = exp£p e (r c>1 + r c>2 ) - + x’ ^ - n T (y' cl + 


(5) 


where r c ^ and r c 2 are the respective distances from the point of pair production on 
the Z-axis to the surface of the center Nal crystal, x’ . and x’ 0 are the respective 

C, X ^ 3 ^ 

thicknesses of wall material traversed, and y' % and y' 0 are the respective thick- 
nesses of reflector material traversed for the annihilation quanta 1 and 2. These thick- 
nesses are dependent on the angle 0 for the emission of the annihilation quanta. This 
dependence is summarized in the following relationships: 


if e<j, 


•X 


Rr 



sin 0 

t 

x c 

C ,1 

sin 0 

T 

y c 

C ,1 

sin 0 


\ for tan 0 = 


R r 


L - t 


(6a) 


j 
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nN 


[ ’c,2 = 


x' 


Cj 2 

y' c ,2 = 


Rc 

R c1 

sin(0 + 7 r) 

sin 9 

x c 

x c 

sin(0 + 77 ) 

sin 9 

y c 

y c 

sin(0 + 77 ) 

sin 6 
J 


R r , R r 

^ ^ — i- / <n\ > C 


> for tan(0 + 77) = -^ or tan(0) = -p 


(6b) 


■c,2 = 


X c,2 


c,2 


t 

t 

COS (0 + 7 r) 

cos 9 

X C 

x c 

cos(0 + 7r) 

cos 0 

y c 

y C 

cos(0 + 7 7) 

cos 0 
J 


Rr 


\ for tan(0 + 77) = tan 0 < — 


(6c) 


R, 




C >1 sin 0 


x 


c,l sin 9 


, = Jc_ 

^c,l sin 0 


R r 


> for tan 0 = -tan 0 > — — 


(6d) 


'c,! 


t 

cos 0 


t 




cos 0 


C c,l 


cos 6\ 


cos 9 


> for tan 9 


= -tan 9 = 


Rr 


(6e) 


'c,l 


y c 

cos 9 


y c 

cos 9 
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r c,2 = 


R c 

_ R c 

sin(0 + tt) 

sin 0 

x c 

- Xc > 

sin(0 + tt) 

sin 0 / 

y c 

y c 

sin(0 + 7 t) 

sin 0 
> 


R 

for |tan(0 + ir)| = -tan 0 = — ^ 


(6f) 


where L is the length and R c is the radius of the center Nal detector. If 8 < ^ and 

R c 2 

tan 6 < , annihilation quantum 1 will escape from the center crystal through the 

L - t 

photomultiplier tube and has no chance of being captured by the side detectors. Simi- 

T> 

larly, if 8 > ^ and |tan(0 + -n) = -tan(0) < ^ annihilation quantum 2 will escape 

being captured by one of the side detectors. In either of these cases, the probability of 
capturing both annihilation quanta with the side detectors is zero. 


Processes in Annihilation-Quanta Detectors 

If both quanta of a pair are able to escape from the center detector and its housing, 
the probability that one, or both, will interact with the side detectors can be predicted by 
considering the direction of emission, the attenuation by the housing and reflector of the 
side detector, and the interaction of the quanta with the Nal crystals. If a quantum (for 
example, quantum 1) enters one of the Nal side detectors, the probability that it will inter- 
act with the crystal is 


Pj k (interaction of l] = 1 - exp 

where r g j is the distance that a quantum traverses through a side detector along a path 
with origin at O'. This distance is given by 


-p e r 


e 1 s,l 


( 7 ) 


r s,l =*1 - 


P 

sin 0 


(8) 


where q^ is the total distance from O' to a point on the surface of the side detector as 
determined by 0 and cp, and p is the radius of the hemicylindrical groove cut into the 
face of the detector. The same relations hold if cp is such that annihilation quantum 1 
passes through the other side detector. Because of this symmetry about the Y-Z plane, 
the remaining discussions will assume that annihilation quantum 1 is emitted only when 
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-?- = cp = and in the final summation over cp the resulting probability will be multi- 

plied by two. Similarly, the probability that the other annihilation quantum will interact 
with the Nal crystal of the other side detector is given by 


Pj ^[Interaction of 2] = 1 - exp 


-M e r s,2 


(9) 


where r g 2 is the distance that annihilation quantum 2 traverses through the crystal 
along a path with origin at O' and is given by 


r s,2 = *2 + 


sin(0 + u) 


q 2 - 


sin 0 


(10) 


where q 2 is the total distance from O' on the Z-axis of the center detector to a point 
on the surface of the side detector as determined by 0 and cp. 

The following discussion of determining the distances q^ and q 2 is more easily 
visualized by redefining the coordinate system for the side detectors as shown in figure 5. 
For simplicity the grooves in the face of the detector are omitted. The £-axis is along 
the Z-axis in figure 4. If 0 and cp are such that annihilation quantum 1 exits through 
the flat face of the upper crystal as in figure 5(a), 


cot 0' = cot 9 + — cos (p 

H r 


(ID 


where d is the distance of the point of annihilation pair production on the £-axis from 
the origin, and H is the crystal height. The position d is given in terms of the orig- 
inal X,Y,Z coordinate system by 


d = (I + t) - R s (12) 

where I is the indentation of the center detector into the side detectors, R s is the 
radius of the side detectors, and t is as defined in figure 4. 

On the face of the crystal, 

n 2 + c 2 = Rg (13) 

where 


rj = sin 0'sin cp 
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£ = Tj cos 0' 


Also, 


or by substitution, 


Ti = 


H 

sin 0'cos cp 


H 2 tan 2 cp + H 2 S Rg 

COS 2 


( 14 ) 


By substituting equation (11) for cot 6 ' into equation (14), the necessary condition for 
the quantum to exit through the face of the side crystal is given by 




cot 0+^7 COS cp 
0 v H 

ia.rr'cp + 


cos 2 cp 


= Re 


(15) 


If this relation holds for a particular value of 0 and cp, then qj is given by 


Ql = 


H 

sin 0 cos cp 


(16) 


Because of the symmetry about the rj - £ plane, the necessary condition for the other 
quantum to exit through the end face of the other crystal is given by 


w2 


t-o n 2 f n 
T 


cot 0 - ■£ cos cp 

COS 2(p 


(17) 


If the condition of equation (17) is met, then 


q 2 - 


H 

sin 0 cos cp 


(18) 


If either condition (16) or (18) is not met, one or both of the annihilation quanta will exit 
through the side of the crystal; and the distances from the point of creation through either 
crystal can be found by using figure 5(b) and the relation 

T 1 = Rg + i 2 (19) 
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where 


and 


ij = q^ sin 8 cos cp 


T 1 


= Ql 


sin 6 
sin O' 


From the triangle formed by the sides d, q^, and r^, it can be shown that 


cot 6' = cot 6 + 

q^ sin 0 


(20a) 


and squaring yields 


cot 2 0 ' = cot 2 9 ♦ M £2^9 + <L_ 

^1 sin^0 q|sin^0 


(20b) 


From equation (19) it can be found that 


R‘ 


csc 2 0' = s + COS 2 (p 

q 2 sin 2 0 


( 21 ) 


By using the trigonometric relationship, csc2 0 ' = 1 + cot 2 #', equations (20) and (21) can 
be combined to yield 


1 + cot 2 9 ♦ f 22^ + = + cos 2 v 

*1] sin' 2 # q"sim0 q"sin^0 


( 22 ) 


or upon rearranging 


q^(l - sin 2 0 cos 2 <^ + 2q^d cos 0 + ( d 2 - R^) = 0 


or solving for q^ 


-d cos 0 + 

d 2 cos^0 + | 

" r 2 ,2' 
[ S J 

) ^1 - sin 2 0 cos 2 <^j 

1 

*1 - sin20 cos2 cpj 


qj = 


The positive radical is used because q^ cannot be negative. 


(23) 
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Again by using the previous symmetry arguments, q2 can be found to be 



Under actual operating conditions the single -channel analyzers connected to the side 
detectors have their lower levels and window widths adjusted to accept only pulses with 
amplitudes vnthin the width of the full-energy peak for the 0.511-MeV photons. This 
requirement makes certain that only pair-production processes are analyzed; that is, 
Compton scattered photons from the center detector probably will not be able to trigger 
the coincidence unit. Therefore, the interaction probability for the side detectors must 
be multiplied by the peak-to-total counting ratio (called photofraction) for the particular 
detectors used. Theoretical data presented in reference 5 indicate that for Nal detectors 
with diameters from 7.6 cm to 12.7 cm and heights from 7.6 cm to 12.7 cm the peak-to- 
total ratios lie within the range 0.75 to 0.80 for the 0.511-MeV photon. By accounting for 
attenuation by the light shield and reflector material of the side detectors and by assuming 
that one annihilation quantum is captured by each of the side detectors, the probability of 
capturing both of the quanta which escape from the center detector within the photopeak of 
the side detectors is given by 


k [pair captur^ = iexp - f x s + p r y s ^ \ F 2 Pj interaction of l)Pj k (Tnteraction of 2] (25) 


where /if and /i r are the respective attenuation coefficients for the light shield mate- 
rial and reflector material of the side detectors, x s is the shield thickness, y s is the 


reflector thickness, and F is the peak-to-total ratio for the side detectors for the 


0.511-MeV photons. 


By combining the probabilities for pair production, escape from the center detector, 
attenuation by the detector housings, and capture of pairs, and by summing over 8 and 
< p , the total probability (or efficiency) of producing and capturing the annihilation quanta 
for an incident photon with energy E is given by 
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The summation over i is terminated at 

t = 2R S - I. (27) 

since for t < 2R S - I, at least one photon of a pair would escape being captured by a side 
detector. The summation over k allows cp to vary from - A cp to - A <p so 

that there will be no contribution for <p = — or for cp = The reason for this limita- 

U Lt 

tion is that the two side-detector crystals are separated by approximately 0.6 cm. With 
the proper choice of a the loss in available solid angle for capturing the pairs may be 
taken into account. 


RESPONSE OF THE SPECTROMETER 

As stated previously, in order to use a spectrometer such as described in this 
report, it is necessary to know the response of the spectrometer to monoenergetic pho- 
tons over the energy range of interest; that is, it is necessary to know the pulse-height 
distributions produced by photons with energy between E and E + AE over the entire 
energy range of interest. If a set of monoenergetic sources existed with energies in each 
increment AE, experimental measurements could be made to furnish the necessary 
response. Since no complete set exists over the energy range from 1.0 MeV to 10.0 MeV, 
the response was calculated by using one available source and the predicted efficiencies 
for the other energies. The simplicity of the pulse-height distribution taken with the 
three-crystal pair spectrometer for a monoenergetic photon source greatly simplifies 
this task. 

As shown in figure 3(b), the pulse-height distribution for the 4.43-MeV photons con- 
sists essentially of the sum of a Gaussian shaped peak and a triangular -shaped low- 
energy tail. In order to calculate the response of the detector to other energies, essen- 
tially the same shapes were assumed; however, the following criteria were used: 

(1) The full widths at half maximum of the Gaussian shaped peaks were made pro- 

portional to the square root of E - 1.022 MeV 

(2) The areas under the triangular -shaped regions were made proportional to 

E - 1.022 MeV 

(3) The total area under each pulse-height distribution was normalized to the total 

efficiency of the detector at the particular photon energy 

The first criterion is based on the usual convention that the resolution of a scintillation 
detector is inversely proportional to the square root of the energy deposited in the 
detector, which for a pair -production process is E - 1.022 MeV. Consequently, the 
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actual width would be proportional to the square root of the energy deposited. The sec- 
ond criterion arises from the use of the approximation that the fraction of electron or 
positron energy converted into bremsstrahlung in a thick target is proportional to the 
incident energy. (See ref. 6.) The constants of proportionality for the Gaussian widths 
and areas under the triangular- shaped tails were determined from the measured 
4.43-MeV-photon pulse-height distribution. This procedure was programed, and a 
63 by 63 response matrix was computed. 

UNFOLDING OF THE PHOTON SPECTRA 

The measured pulse-height distribution taken with a spectrometer, as described in 
this report, of a multienergetic photon source or a continuous distribution of photon 
energies is a result of summing the pulse-height distributions produced by a number of 
monoenergetic photons. Because of the finite resolution o^the spectrometer and the 
bremsstrahlung and electron losses which are dependent on the incident photon energy, 
the instrument response is nonunique; and some method of analyzing, or unfolding, the 
pulse-height distributions must be used to give the actual incident photon spectra. 

If an incident photon spectrum is represented by a column matrix N and if the 
response of a spectrometer is represented by a matrix B, the pulse-height distribution 
M produced by the spectrometer is given by 

M = BN (28) 


The problem is then reduced to solving the matrix equation 

N = B" 1 M (29) 

This straightforward approach has been found to yield poor results because inherent 
errors in the response matrix are magnified by the inversion process. (See ref. 7.) 
Therefore, the method used in the present report is the iterative method described in 
reference 4, which is discussed in some detail in the appendix. This iterative procedure 
has been shown to yield much better accuracy than the matrix inversion method. (See 
ref. 7.) 


RESULTS AND DISCUSSION 

The dimensions of the spectrometer are given in table I. The attenuation coeffi- 
cients used for determining the spectrometer efficiency were obtained from reference 8 
and are given in table II. 
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Equation (26) for the total probability of detection (or efficiency) has been pro- 
gramed for the digital computers at Langley Research Center, and the calculated efficien- 
cies are given in figure 6 by the solid line drawn over the photon energy range from 
1.022 MeV to 10.0 MeV. Below 1.5 MeV, the pair -production cross section for Nal is 
very small; therefore, experimental verification of the efficiency of the detectors at these 
low energies would require very long counting times and be subject to relatively large 
statistical errors. 

The only readily available source of relatively high-energy (greater than 3 MeV) 
gamma rays was an americium-beryllium (Am-Be) neutron source which gives off a 
4.43-MeV gamma ray. The activity of this source was calibrated by measuring the num- 
ber of 4.43-MeV gamma rays in the photo peak of spectra taken with a 7.62-cm by 
7.62-cm Nal scintillation detector. The efficiency and photofraction for the detector were 
obtained from reference 9. By using this source, the measured efficiency of the three- 
crystal pair spectrometer was jpund to be 7.9 x 10"^ as shown by the circular symbol in 
figure 6. The error bars indicate the combined possible error due to counting statistics 
for both the single-crystal detector and pair spectrometer and placing an uncertainty of 
±10 percent on the data of reference 9. The square symbols shown in figure 6 are effi- 
ciency calculations found by using only simple pair-production attenuation, and neglecting 
geometrical considerations. Equation (26) then reduces to 



where fi g is the pair-production attenuation coefficient at photon energy E, R s is the 
radius of the side detectors, I is the indentation of the center crystal into the well, and 
e is a normalization factor. The normalization factor was obtained by equating the cal- 
culated efficiency of the system from equation (26) at 5 MeV and equation (30). At 10 MeV 
the simple calculation differs by 13 percent from the calculations which include the spe- 
cific geometry of the detector system. If a normalization point had been chosen with 
lower photon energy, the relative difference would be increased for the higher photon 
energies. If the spectrometer efficiency were experimentally measured with a calibrated 
source such as the Am-Be source with the 4.43-MeV gamma ray, the use of the simple 
calculation (eq. (30)) to obtain the efficiency at higher energies would cause an underesti- 
mation of at least 13 percent at 10 MeV. 

After a Compton scattering, the scattered photon can possibly undergo a pair- 
production process. Calculations have been made to determine the error introduced by 
neglecting this type of multiple interaction. These calculations were made by finding the 
average energy of Compton scattered photons produced by 10.2-MeV and 6.6-MeV incident 
photons on Nal and by finding the average depth in the center detector of such Compton 
scatterings. The probability of having such an average scattering and then having the 
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scattered photon undergo a pair-production process in the remaining detector thickness 
was then computed. This probability was found to be 4.5 percent and 3.8 percent of the 
probability of having single pair-production processes at incident energies of 10.2 MeV 
and 6.6 MeV, respectively. These percentages should be an overestimation because no 
allowance was made for the Compton scattered photons being deflected from the incident 
direction. 

The stripping approximation and the iterative unfolding procedure were programed 
for the Control Data series 6000 computer system at Langley Research Center, and a test 
case was run to check the convergence of the iterative procedure. A simulated brems- 
strahlung spectrum was multiplied by the response matrix to yield a simulated pulse- 
height distribution. This simulated pulse-height distribution was stripped to yield the 
first approximation to the energy spectrum and then unfolded by the iterative procedure 
to yield the approximate bremsstrahlung energy spectrum. Table III gives the input 
energy spectrum, the simulated pulse -height spectrum, the stripped approximation, the 
first approximation as used in reference 4, and the iterative approximation to the energy 
spectrum after 201 iterations. As shown, the iterative approximation is within 1 percent 
except for photon energies of 10.35 MeV and 10.50 MeV. It can also be seen that the 
stripping procedure approximates the input much better than does the first approximation 
of reference 4. This improvement should greatly reduce the number of iterations neces- 
sary for good convergence. 

The improvement obtained by using the iterative procedure rather than inverting 
the response matrix as in equation (29) can be seen immediately in figure 7, which shows 
the percent difference between the input spectra of table III and the unfolded spectra by 
both the iterative method and the matrix inversion method. The data points for the 
matrix inversion method are not shown for channels at which they did not differ from the 
iterative method sufficiently enough to be distinguished in the figure. The matrix inver- 
sion process is shown to exhibit oscillations about the input data, as mentioned in refer- 
ence 7. These oscillations or instabilities generate errors as large as 30 percent. 
Although these oscillations for continuous spectra may be smoothed out, the direct inver- 
sion process is not well suited for spectra containing discrete energy peaks. 

CONCLUDING REMARKS 

An investigation has been made of the use of a three-crystal spectrometer for mea- 
surements of photon spectra. Calculations have been made of the efficiency and response 
of the spectrometer system for a photon energy range of 1.5 MeV to 10.0 MeV. An itera- 
tive spectrum -unfolding method has been modified to permit operation with a spectrom- 
eter system of low efficiency, such as the three-crystal spectrometer. The efficiency 
calculations appear to have good accuracy when compared with a limited amount of 
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experimental data and would appear to be an improvement on simple calculations which 
neglect the specific geometry of the system. The iterative unfolding technique has proved 
to yield errors which are insignificant in comparison with errors in either efficiency or 
response and has proved to be much superior to unfolding by the matrix inversion method. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 19, 1969, 

124-09-11-03-23. 
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APPENDIX 


ITERATIVE ANALYSIS OF PHOTON SPECTRA 

The method to be described and used to unfold the photon spectra is based on the 
iterative method described by Smith and Scofield in reference 4. The incident photon- 
energy spectrum is represented by a column matrix N in which each row element N n 
represents the number of photons in an energy interval n. The measured pulse-height 
distribution is represented by a column matrix M in which each row element M m 
represents the number of pulses in a pulse-height interval m. The response of the 
detector to some photon with energy in the interval n can be represented by a rectangu- 
lar matrix B in which the elements B m n of the nth column are the probabilities that a 
photon with energy in the nth energy interval will produce a count in the mth pulse-height 
interval. The problem of analyzing some measured pulse-height distribution is then 
reduced to finding the elements N n where 

M = BN (Al) 

If equation (Al) is multiplied by the transpose of the response matrix, then 

BM = BBN or S = C N (A2) 

where the matrices S and are defined as 

S = BM and C = B B 

The matrix C is square and symmetric. 

The iterative method of reference 4 consists of calculating a diagonal matrix D 
as follows: 

N = D S (A3) 

Since D is diagonal, S and N must be of the same order, and equation (A3) can be 
written 


Np - DppSp (A4) 

where Np are the elements of N, Sp are the elements of S, and Dpp are the diag- 
onal elements of D. 
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APPENDIX — Continued 
The elements of D can be found by 

Nn , , 

D pp - (A5) 

b P 


where the off-diagonal elements, of course, are zero. The iterative scheme consists of 
the following steps: 

(1) A set of values for the energy spectrum N q is approximated, where q 
represents the number of iterations, and q = 1 for steps (1) through (4) 

(2) Approximate values for the S matrix are calculated by 


S q = C N q 


(q = i) 


(A 6) 


(3) An approximate set of diagonal elements for D is computed by 


D q = — 2 

pp s q 

p 


(q= 1) 


(A7) 


(4) The sum of the squares of the differences of the elements of S = BM and S q 

are computed by 

T q -(s p -s^ 2 (q = 1) (A8) 

(5) The next approximation for N q is computed by 

N q+1 = D q S (A9) 


(6) The next approximation for 


(7) The next approximation for 


S is computed by 
S q+1 = C N q+1 


D is computed by 


, N q+1 
D q+1 = P 
PP .q+l 
S P 


(A10) 


(All) 
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APPENDIX — Continued 

(8) The test for the convergence of is made by 



(9) Steps (5) to (8) are repeated until 

t A- 1 _ T q < T q - T q+1 (A13) 

or until q reaches a predetermined value 

The computed when either of the conditions of (9) is met is considered to be the 
incident photon-energy spectrum. A thorough discussion of the convergence and unique- 
ness of the solution is given in reference 7. 

The remaining problem to be solved is to make a reasonable first approximation for 
N. In reference 4 this approximation was accomplished for the case of a single-crystal 
Nal detector by letting 

N 1 = S = B M (A14) 

However, the detection efficiency of the three-crystal pair spectrometer is several 
orders of magnitude less than that of a single-crystal Nal spectrometer; therefore, the 
use of equation (A14) would result in a severe underestimation for the first approximation 
of the iterative procedure. Because of this underestimation a somewhat different 
approach is used herein. 

A method of stripping the measured pulse-height distribution into a photon-energy 
spectrum was used for the first approximation for N for the three-crystal spectrom- 
eter. Because of the narrow Gaussian peaks in the response matrix, the last element of 
the pulse-height distribution M u is mainly a result of the last two elements of the pho- 
ton spectra N v and N v _j. Consequently, 

M u ~ ®u,vNv + Bu^y-lNy-l (A15) 

If N v is approximately equal to N v _j, then 



B u,v 


M 

B u,v 


1 


(A 16) 


The pulse-height distribution can be stripped by subtracting the contribution to each 
element of the pulse-height distribution by the number of incident photons with energy in 
the vth energy interval from the elements of M: 
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APPENDIX - Concluded 
M m ~ M m - B m>v N v 


Similarly, 


N v -i 


M. 


v - 1 


*u-l,v-l + B u-1,v-2 


By continuing these steps, a general recursion relation can be found so that the 
photons in the nth energy interval can be approximated by 

u 

^n - | B n ,mNm 

m=n+l 

n ®n,n + ®n,n-l 

where n = v, v-1, v-2, ... 1 and equation (A 16) holds if n = v. 


(A 17) 


(A 18) 
number of 


(A 19) 
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TABLE I.- SPECTROMETER DIMENSIONS 


Center crystal: 

Length 10. 16 cm 

Radius 2.06 cm 

Light-shield thickness 0.04 cm 

Reflector thickness 0.28 cm 

Side crystals: 

Height 7.62 cm 

Radius 6.35 cm 

Light-shield thickness 0.05 cm 

Reflector thickness 0.11 cm 

Radius of groove 2.70 cm 

Indentation of center detector into well 2.86 cm 


TABLE H.- ATTENUATION COEFFICIENTS 


Photon 

energy, 

MeV 

Attenuation coefficients, cm“l 

Nal 

total 

Nal 

pair production 

A1 

total 

Fe 

total 

A1 2°3 

total 

0.5 

0.350 

0.0000 

0.228 

0.659 

0.300 

1.0 

.215 

.0000 

.166 

.469 

.219 

1.5 

.172 

.0027 

.135 

.381 

.178 

2.0 

.151 

.0090 

.117 

.333 

.152 

3.0 

.134 

.0234 

.095 

.283 

.126 

4.0 

.129 

.0369 

.084 

.259 

.109 

5.0 

.127 

.0482 

.076 

.246 

.098 

6.0 

.128 

.0577 

.072 

.239 

.088 

8.0 

.130 

.0731 

.065 

.233 

.082 

10.0 

.135 

.0894 

.062 

.232 

.077 

15.0 

.147 

.1115 

.059 

.240 

.070 
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TABLE EL- UNFOLDING- PROCEDURE TEST CASE 



30 000 
28 000 










29 


igure 1.- Three-crystal pair spectrometer. 







PAIR SPECTROMETER 



Figure 2 .- Spectrometer electronics diagram. 























Figure 6.- Efficiency of three-crystal pair spectrometer. 
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Figure 7.- Difference between input spectrum and unfolded spectra. 
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